A repulsive trap for two electrons in a magnetic field 
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We study numerically and analytically the dynamics of two classical electrons with Coulomb in- 
teraction in a two dimensional antidot superlattice potential in the presence of crossed electric and 
magnetic fields. It is found that near one antidot the electron pair can be trapped for a long time 
and the escape rate from such a trap is proportional to the square of a weak electric field. This 
is qualitatively different from the case of noninteracting electrons which are trapped forever by the 
antidot. For the pair propagation in the antidot superlattice we found a broad parameter regime 
for which the pair is stable and where two repulsive electrons propagate together on an enormously 
large distance. 
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I. INTRODUCTION 

Recent technological development allowed to create 
various types of surface superlattices for two-dimensional 
(2D) electron gas in semiconductor heterostructures and 
to investigate their transport properties in presence of 
magnetic field. Experiments with antidot lattices VffP 
carried out by different experimental groups (see e.g.Efu) 
and the contribution of classical periodic orbits in the 
resistivity peaks at certain values of magnetic field was 
clearly identified. In fact the size of antidots and the 
distance between them are relatively large, and analysis 
of classical trajectory dynamics can be successfully ap- 
plied to understand a number, of unusual transport prop- 
erties in such antidot arraysB'tl Due to nonlinearity of 
motion in the vicinity of antidot potential the classical 
dynamics can be chaotic that leads to diffusive spr< 
ing of trajectories even for perfectly periodic lattice, 
If the distance between antidots is large or comparable 
to the cyclotron radius of an electron in a magnetic field 
perpendicular to the lattice, then one should understand 
first, the properties of electron dynamics near one an- 
tidot. In the absence of electric field, the dynamics is 
integrable for an antidot of circular shape due to angu- 
lar momentum conservation and an electron always regu- 
larly rotates around the antidot. An electric field applied 
in the 2D plane of the superlattice breaks the cylindri- 
cal symmetry and can lead to the electron escape from 
the antidot to infinity. The problem of electron dynam- 
ics in crossed electric E and magnetic B fields near a 
circular elastic disk (antidot) was studied intu. It was 
shown that the dynamics can be described by a simple 
area-preserving map which depends only on one dimen- 
sionless parameter e = (2Trm/ae)(E / B 2 ), where m, e are 
electron mass and charge, and a is the disk radius. For 
small e < e c the electron dynamics in the phase space of 
angular momentum I and conjugated angle 4> is bounded 
by the invariant Kolmogorv-Arnold-Moser (KAM) curves 
so that electron always remains near the disk. On the 
contrary for e > e c the KAM curves are destroyed, global 



chaos sets in and the electron escapes to infinity after few 
collisions with the disk. 

Until now the classical dynamics in antidot lattices was 
studied only for noninteracting electronsBiJ. In this pa- 
per for the first time we analyze the effect of Coulomb 
interaction between classical electrons in the vicinity of 
an antidot. We show that for sufficiently strong interac- 
tion between electrons their dynamics becomes chaotic. 
Due to that one or two electrons can escape from the 
antidot even in an arbitrary weak applied electric field 
E that corresponds to e c — > contrary to e c > in ab- 
sence of interaction. We determine the dependence of 
average escape rate Y on e showing that in the limit of 
small electric field Y oc e 2 . After that we also discuss the 
two electron propagation in the antidot superlattice. 

The paper is organized as follows. In the next Section 
we briefly discuss the one electron dynamics near the an- 
tidot in crossed magnetic and electric fields. In Section 
III the dynamics of two interacting electrons is analyzed 
in detail. The electron motion in the antidot superlat- 
tice is considered in Section IV. In the last Section we 
summarize the obtained results. 



II. ONE ELECTRON DYNAMICS 

The dynamics of an electron in crossed electric and 
magnetic fields in two dimensions with one antidot is 
described by the Hamiltonian: 



H = (p - eA) 2 /2m + U(x, y) - eEr 



(1) 



where A is the vector potential and U (x, y) describes 
the antidot potential which depends only on the radius 
r = \Jx 2 + y 2 with U = for r > a. For conve- 
nience, followingEl, we introduce the dimensionless vari- 
ables x = x/a,y = y/a,t = t/ro,Ho = Ho/2eF,U — 
U/2e F , B = B/B , E = E/E where e F (v F ) is the Fermi 
energy (velocity), To = (2ej? /ma 2 ) -1 ' 2 = cl/vf, magnetic 
and electric fields are scaled by Bq = (mtp) 1 ! 2 / 'ea and 
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Eq = 2eF / ea respectively. In these units a magnetic field 
B = Bo gives the cyclotron radius R c — a for electron 
with energy ep = mvp/2. We choose the Landau gauge 
A = (—By, 0,0). Then omitting the tildes the Hamilto- 
nian equations of motions reads: 

dx/dt — v x , dvx/dt — Bv y — dU/dx — E x , (0 \ 
dy/dt — Vy, dvy/dt = —Bv x — dU/dy — E y 

where v x = p x + yB, v y = p y . To model the antidot we 
chose the potential: 

U(x iy ) = U (l-r) 6 . (3) 

Usually we choose Uq to be much larger than the electron 
energy H so that this potential becomes very similar to 
an absolutely rigid disk with effective radius a e f f about 
15% smaller than a. 

Far from the antidot the equations of motion are ex- 
actly solvable and give electron rotation over a circle of 
cyclotron radius R c — v/uj c with cyclotron frequency 
lu c = B. In addition this circle moves with the drift ve- 
locity Vd = E / B in the direction perpendicular to electric 
and magnetic fields. As it was found by Berglund et al.u, 
near the antidot, the dynamics strongly depends on the 
dimensionless parameter e = v c {2,i:/uj c = 2irE/B 2 . For 
e«l the electron scatters on the antidot and escapes to 
the infinity after one collision. On the contrary, the sit- 
uation with not very large e is much richeru. In this case 
electron can collide many times with the antidot and this 
process is described by a simple area-preserving mapD: 

4> = 4> + n — 2 sin -1 /3 

/3 = /3-esin0 [ ' 

where bars denote the new values of variables after colli- 
sion, 4> is the scattering angle measured in respect to the 
direction of drift velocity and is the scattering impact 
parameter divided by the antidot radius. In this way 
(3 varies in the interval (-1, 1). We note that (3 can be 
also considered as the orbital momentum I of the elec- 
tron divided by the maximal momentum l max = o,v at 
which electron still collides with the antidot. The real 
dynamics is correctly described by the map if R c 3> 1, 
that corresponds to u > 5. For e ^> 1 the variation of 
/3 is bounded by the invariant KAM curves and electron 
is trapped near the antidat. The last KAM curve is de- 
stroyed for e > e c « 0.45 13 so that the orbits with initial 
[3 w can escape from the antidot to infinity. Of course, 
for e > e c some islands with regular motion inside still 
remain, but they become very small as soon as e becomes 
significantly lager than e c . 

To study the electron dynamics in ([!]) the Hamiltonian 
equations of motion are solved numerically by Runge- 
Kutta method of fourth order so that the electron energy 
is conserved with the relative precision better than 10~ 6 . 
The examples of the Poincare cross sections constructed 
at x = and v x > for trajectories trapped near the 
antidot is shown in Figs. 1,2. In Fig. lew 0.16 is rather 



small and almost all phase space is filled by integrable 
KAM curves. For Fig. 2, the parameter e « 0.42 is 
close to e c and KAM curves become more deformed and 
the chaotic component becomes visible. This case can be 
compared with Fig. 3 inEI where the cross section for the 
map (^) is given for a close value of e. 



Py 




Y 



FIG. 1. Poincare cross section for the Hamiltonian (|l|) con- 
structed at x = 0, v x > for H = 8.725, B — —2, E = 0.1 so 
that e ~ 0.16. The antidot determined by the potential (^) is 
located at (0, 0); U = 1000. 
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FIG. 2. Same as in Fig. 1 but for E = 0.3, e w 0.42. 

It is interesting to note that the map (Q) can be not 
valid if orbits have \j3\ w 1 or |/3| > 1. For example, if 
the antidot is inside a large cyclotron circle then elec- 
tron will make many many rotations before this slowly 
drifting circle will cross the antidot. This situation is not 
taken into account by the first equation in (Q) . An exam- 
ple of the electron dynamics in this case is given in Fig. 
||. Here e ~ 0.16 is small and the motion is still regular. 
We should stress that such type of trajectories separates 
orbits which escape to infinity and those which collide 
with the antidot on each cyclotron period. For the an- 
tidot superlattice with antidot spacing comparable with 
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R c this type of orbits (see Fig. ||) is of special importance 
since these orbits can easily jump from one antidot to an- 
other leading to a global diffusion in the system. We will 
discuss this situation in Section III. 
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FIG. 3. Electron dynamics in (x, y) plane near the antidot 
of Fig. 1 for H = 9.7, B = -2.0, E = 0.1, e « 0.16. 



III. EFFECTS OF COULOMB INTERACTION 
ON ELECTRON DYNAMICS 

Let us now consider how the Coulomb interaction be- 
tween two electrons affects their dynamics near the anti- 
dot. In this case the Hamiltonian of the system reads 

H = i/ (Pi,ri)+ff (P2,r 2 ) + e 2 /|ri-r 2 | . ( 5 ) 

While in the free space the Coulomb interaction repels 
the electrons and leads to their separation the situation 
is more complicated in the presence of magnetic field. In 
the case without any antidot the total momentum of two 
electrons is conserved and as a result each electron rotates 
regularly on a cyclotron circle which in addition rotates 
around the center of mass of the system. Without an ex- 
ternal electric field (E — 0) the center of mass is fixed and 
inert whereas in the presence of the field (\E\ > 0) the 
center of mass drifts with the constant velocity Vd = E/B 
but the average distance between electrons remains con- 
stant. However this electron pair can be trapped by the 
repulsive potential of the antidot so that the electrons 
will spend a long time colliding with this antidot. An ex- 
ample of electron dynamics in this case is shown in Fig. 
4. It shows that an electron can escape from the antidot 
even in the situation with e < e c when without the in- 
teraction the electrons remain trapped near the antidot. 
In our numerical simulations we observed different cases 
where one or both of the electrons escape to infinity. 
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FIG. 4. Dynamics of two electrons in {x, y) plane for 
H = 15.65, B = -2.0, E = 0.15, e « 0.24 and initial distance 
between electrons |ri — r2| ~ 0.5. After many cyclotron pe- 
riods the first electron escapes from the antidot to infinity 
(upper figure) while the second remains trapped forever (bot- 
tom figure). 

To investigate how the escape rate depends on the 
strength of an external electric field E we studied the en- 
semble of 100 paths. In each path the positions and mo- 
mentums of each electron are chosen randomly in the in- 
tervals —4 < x, y < 4, —2 < p Xl p y < 2 in such a way that 
the total energy is H w 15±0.5. We remind that the an- 
tidot with potential (j|) is placed at (0, 0) and U$ = 1000. 
The escape rate T is defined as V = 1/T where T is the 
time after which the distance of one of the electrons from 
the antidot is greater than R esc « 5R C 10. This dis- 
tance is sufficiently large and as soon as it is reached an 
electron escapes to infinity and never returns back to the 
antidot. The average value of T is obtained by averaging 
over 100 values obtained for 100 randomly chosen paths. 

The dependence of the escape rate T on the strength 
of the applied electric field E is presented in Fig. 5. It 
definitely shows that the escape takes place even at very 
weak electric fields with e <C e c when without Coulomb 
interaction electrons are forever trapped near the antidot. 
According to the obtained numerical data (see Fig. 5) in 
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the limit of e — > the escape rate is 



(0) 



Our understanding of this dependence is based on the 
following argument. Due to the Coulomb interaction be- 
tween electrons their dynamics in the vicinity of the an- 
tidot becomes chaotic. Therefore the phase (j> in the map 
([|) changes randomly between electron collisions with 
the antidot and (3 grows diffusively with the number of 



collisions n so that (A(3) ~ Dn with D = e 2 /2. 
diffusion results in the escape rate T/uj c ~ D ~ e 2 
in agreement with the numerical data in Fig. 5. 
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FIG. 5. Dependence of the escape rate V on the electric 
field E for electrons being initially in the vicinity of the anti- 
dot at B = -2.0, H ps 15.0 ± 0.5, U = 1000. Here averaging 
is done over 100 paths, uj c = 2, e = 2ivE/B 2 , points give the 
numerical results for V and the straight line gives the depen- 
dence (^|). Logarithms are decimal. 



IV. TWO ELECTRON PROPAGATION IN 
ANTIDOT SUPERLATTICE 

We also studied the electron dynamics on a square an- 
tidot superlattice, when the antidot potential is given by 
(|J) and the distance between antidots d > 2. In this case 
our results for one electron dynamics are. in qualitative 
agreement with the conclusions drawn ind'EI. As soon as 
the cyclotron radius R c becomes comparable with the an- 
tidot spacing d the trajectories start to move diffusively 
on the whole lattice. Trapped orbits near one antidot 
exist only for 2R C < d and e < e c . For R c > d/2 the 
cyclotron circle starts to drift in a way similar to that 
one shown in Fig. 3. After a time td ~ d/va a collision 
with another antidot takes place that finally originates a 
sequence of irregular jumps between antidots. The dif- 
fusion rate on the superlattice originated by this process 
can be estimated as Di at ~ d 2 /td ~ d Vd- This diffusion 
is important in the limit R c ~ d 3> 1. However, we note 
that even at E — at R c > d/2 there are chaotic orbits 



which diijjuse over the whole lattice as it was discussed in 
detail inQij and this diffusion is dominant for d ~ 1. 

It is interesting to understand how two electrons move 
in such a superlattice. Intuitively, one would expect that 
the Coulomb repulsion will separate electrons and they 
will not propagate together. In fact we found that it is 
not necessarily the case and there are regimes where two 
electrons propagate together. An example of such a case 
is shown in Fig. 6. In this case the electron pair moves 
with an average drift velocity Vd ~ E/B and the total 
displacement of the pair is about hundred times larger 
than the distance between the two electrons (Fig. 6). 
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FIG. 6. Dynamics of two electrons in the (a;, y) plane of an 
antidot superlattice for H « 51, B = -2.0, E = 0.4, e w 0.63 
and initial distance between electrons |n — ra| ~ 4. The 
antidots are placed on a square superlattice with period 4 and 
Uo = 1 in The upper figure shows the propagation of the 
first electron in the plane [x, y), while the bottom figure shows 
the distance between the electrons Aa; = xi— x%, Ay = y\—yi. 

The physical reason for appearance of such electron 
pairs is quite clear in the absence of superlattice poten- 
tial. In this case, as it was discussed before at the be- 
ginning of Section III, electrons rotate one around an- 
other and propagate together, their dynamics is inte- 
grablc. Then according to the KAM theorem a weak 
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perturbation will not destroy such pairs. Indeed, in the 
case of Fig. 6 the antidot potential is relatively weak 
(Uo = 1 <C H/2 w 25) and the pair is not destroyed. We 
checked numerically that provided Uo ~ H/2 the pair 
size starts to grow diffusively due to random scattering 
on a strong antidot potential, and eventually the pair is 
destroyed and electrons continue to propagate separately. 
For Uq ^> H the separation happens after a few collisions 
with antidots. On the contrary for Uo *C H the life time 
of the classical pair becomes infinite in agreement with 
the KAM theorem. 
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V. CONCLUSIONS 



In this paper we investigated the effects of Coulomb in- 
teraction between two electrons on their classical dynam- 
ics in antidot superlattice in crossed electric and magnetic 
fields. We found that for weak electric field the electron 
pair can be trapped for a long time near an antidot even 
if eventually one or two electrons escape from the anti- 
dot. The escape rate F decreases proportionally to the 
square of electric field. This behaviour is qualitatively 
different from the case of noninteracting electrons which 
are trapped forever near the antidot in the limit of small 
electric field. 

The study of the electron pair dynamics in the anti- 
dot superlattice showed that the Coulomb repulsion can 
create stable pairs propagating on a large distance. In 
agreement with the KAM theorem such pairs are stable 
when the antidot potential strength is relatively weak 
compared to the electron energy. On the contrary, in the 
opposite limit the pairs become unstable and electrons 
are quickly separated from one another. On the basis 
of this phenomenon it is possible to make a conjecture 
that in two dimensional heterostructures with high mo- 
bility the impurity potential is relatively weak and such 
electron KAM pairs will be stable and can be detected 
experimentally. We note that in the recent experiments^ 
with 2D electron gas the carriers of charge 2e have been 
detected. It is possible that these carriers are related to 
the KAM pairs found in this paper. 

We thank G.Casati who attracted our interest to the 
results found inQ. 
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